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ABSTRACT 

We report on observations of the Cosmic Microwave Background (cmb) obtained during the January 
2003 flight of Boomerang. These results are derived from 195 hours of observation with four 145 
GHz Polarization Sensitive Bolometer (PSB) pairs, identical in design to the four 143 GHz Planck 
HFI polarized pixels. The data include 75 hours of observations distributed over 1.84% of the sky 
with an additional 120 hours concentrated on the central portion of the field, itself representing 0.22% 
of the full sky. From these data we derive an estimate of the angular power spectrum of temperature 
fluctuations of the cmb in 24 bands over the multipole range 50 < I < 1500. A series of features, 
consistent with those expected from acoustic oscillations in the primordial photon-baryon fluid, are 
clearly evident in the power spectrum, as is the exponential damping of power on scales smaller than 
the photon mean free path at the epoch of last scattering {I > 900). As a consistency check, the 
collaboration has performed two fully independent analyses of the time ordered data, which are found 
to be in excellent agreement. 

Subject headings: Cosmic Microwave Background, Cosmology, Bolometers 



1. INTRODUCTION 

The wealth of cosmological information that is en- 
coded in the statistical properties of the Cosmic Mi- 
crowave Background Radiation (cmb) has motivated 
a highly successful observational effort to measure 
the angular power spectrum of the cmb tempera- 
ture anisotropies. The experimental effort is broad 
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based, with teams reporting results from interferomet- 
ric and single dish observations spanning a decade 
in frequency and more than th ree decades in angula r 
scale [for some recent results, see lBennett et alJ (|2003a|) : 
iBenoit et all i200l : iRuhl et all J20Q31 : lHalverson 

H2004D: II "' 



ll2002tk iReadhead et all il200 J): iKuo et alJ (|2004,. _ 
iDickinson et all l|2004D : iLee et al] l)200lH ]. The success 
of these observations, coupled with the predictive power 
of accurate theoretical modeling, has contributed to the 
ongoing transformation of cosmology into a quantitative, 
precise, and increasingly accurate science. 

In this paper, we present the angular power spectrum 
of the temperature anisotropies derived from the data 
obtained during the January 2003 flight of BOOMERANG 
(hereafter B03). Havi ng been upgrad e d wit h a polariza- 
tion sensitive receiver |Montrov etall l)2003() : iMasi et alJ 
l)2005[) ]. B03 is optimized to probe the polarization of 
the CMB at sub-degree angular scales while retaining full 
sensitivity to the unpolarizcd emission. In addition to a 
measurem ent of the curl-f r ee com ponent of the CMB po- 
larization |Montrov et all ll2005jl| and the temp erature- 
polarization cross correlation |Piacentini et all J2£ 



the 2003 flight has resulted in a precise measurement 
of the angular power spectrum of the temperature 
anisotropics. 

These data represent an improvement on published 
measurements of the temperature power spectrum at 
multipoles 600 < I < 1200, which probe physical scales 
corresponding to the photon mean free path during the 
epoch of last scattering. 

2. INSTRUMENT AND OBSERVATIONS 

Follo wing the successf ul 199 8 Antarctic cam- 
paign |de Bernardis et all l)2000|) : iNetterfield et alJ 
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ll20f)l : iRuhl et all J2001 : ICrill et alJ an en- 

tirely new focal plane was designed around a set of four 
145 GHz Polarization Sensitive Bolometer (PSB) pairs. 
These receivers, the first of their kind to be used in 
astronomical observations, combine the raw sensitivity 
of cryogenic bolometers with intrinsic sensitivity to 
linear polarization, a prope rty his t orical l y associated 
only w i th coherent det ectors |Jonesl l)2005|h Pones et all 

In addition to the four 145 GHz PSB pixels, the 
B03 focal plane accommodates four dual-frequency pho- 
tometers, operating in bands centered near 245 and 345 
GHz. While the CMB temperature anisotropics are de- 
tecte d with high signal to noise in the 2 45 and 345 GHz 
data |Montrovll)2002|klMasi et alJ l)2005|) ]. these data are 
not included in the present analysis. 

2.1. Polarization Sensitive 
Bolometers 

Each bolometer within a PSB pair is sensitive to a lin- 
ear combination of the Stokes I, Q, and U parameters on 
the sky multiplied by the frequency response of the re- 
ceiver, F v , and convolved with the (two-dimensional) co- 
and cross-polar beam patterns, P\\(f) and Pj_(f), respec- 
tively. After deconvolving the system transfer function 
from the time ordered data, each sample from a single 
detector within a PSB pair, di, may be written as the 
sum of a signal component 

* = § / dv w b F» Ifdf (P\\ (r<) + P±(n)) 
I + jV{h) fa cos 2^ + 17 sin ) , 

and a noise contribution. The Stokes parameters are 
understood to be defined on the full sky, and the inte- 
gration variable is f,; = hi — r, for a vector, hi, describing 
the pointing at a time sample, i. The normalized beam 
response and the polarization efficiency are given by, 

m = 7=£§- (2) 

The polarization leakage parameter, e, is defined as the 
ratio of the square of the diagonal elements of the Jones 
matrix describing an imperfect polarizer. For B03's 
PSBs, the leakage ranges from 5% to 8%, as determined 
from pre- flight measurements. The angle ip is the pro- 
jection of the axis of sensitivity of a given detector on 
the sky. This angle is modulated in time due to sky ro- 
tation, and is further affected by the motion of the gon- 
dola. The calibration factor, s, that converts the bright- 
ness fluctuations in I, Q, and U to a signal voltage, is 
obtained through cross-calibratio n with the temperatur e 
anisotropic s obser ved by WMAP |Bennett et"afl <|2003bfl ; 
IMasi et al.lll2005T) 1. 

A mo re detailed d i scussi on of the PSB design may be 
found in lJones et alJ l|2003|) . A description of the charac- 
terizat ion of the B03 PSBs may be found in lMasi et alJ 
l)2005|) . A discussion regarding prop erties and methods 
of a nalysis of PS Bs may be found in iJones et alJ l|20f)5f ) 
and lJoneslpOOl . 

2.2. Sky Coverage 

The Boomerang gondola scans in azimuth, mapping 
the sky signal to a bandwidth between ~50 mHz and 5 
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a The effective beam is defined as the convolution 
of the physical beam with the ~ 2.4' r.m.s. error in 
the pointing reconstruction as determined from point 
source observations in the CMB field. b The focal plane 
noise equivalent temperature, derived from the in- 
flight noise measured at 1 Hz. Note, however, that 
at both low (< 100 mHz) and high (> 5 Hz) frequen- 
cies, the noise rises significantly. The approximate 
noise (~ (diag(Cjv))) per 3.4' pixel in the deep field 
for the data included in this analysis. 

Hz. Sky rotation modulates the orientation of the instru- 
ment with respect to the signal. A complete characteri- 
zation of the Stokes I, Q and U parameters is achieved 
via a combination of this sky rotation as well as the joint 
analysis of the data from all eight detectors in the course 
of pixelization. The low frequency stability of the bolo- 
metric system and balloon environment allow an accu- 
rate reconstruction of all three linear Stokes parameters 
at each pixel with sensitivity to angular scales ranging 
from ~ 10° to the ~ 10' beamsize. 

The launch date, latitude, and the azimuth/elevation 
constraints imposed by solar, limb, and balloon avoid- 
ance considerations limit the region of sky accessible to 
Boomerang. During the Austral summer, the anti- 
solar meridian falls in the vicinity of RA ~ 70° (J2000 
epoch). Continuum emission from the Galaxy is signif- 
icant over much of the available sky; the CMB field is 
chosen to minimize this Galactic contamination, subject 
to the constraints imposed by the LDB flight parameters. 

The Boomerang telescope is an Alt-Az mount that 
is scanned in azimuth at constant elevation. The eleva- 
tion is adjusted on hour time scales. The pendulation 
frequencies of the gondola and scan rate limitations con- 
strain the minimum peak to peak scan amplitude to ap- 
proximately 15°. Earth's rotation provides about 20° of 
cross linking between scans, which aides in the decorrela- 
tion of the Stokes parameters. When the local hour angle 
(modulo twelve hours) is in the vicinity of the RA of the 
target field (i.e., when the sky rotation, and therefore 
the cross-linking, is minimal), our scans are redirected 
to the Galactic plane. 

Given B03's instantaneous sensitivity, the amplitude 
of the CMB polarization signal motivates us to concen- 
trate as much integration time as possible in as small 
a region as possible. However, the competing desire to 
make high fidelity measurements of (the relatively large 
amplitude) unpolarizcd temperature fluctuations led us 
to a two-tiered scan strategy, representing a compromise 
between sensitivity to the temperature and polarization 
power spectra. 

During the first four days of the flight, seventy- five 
hours were dedicated to observations covering ~ 1.84% 
of the sky (these observations are referred to as the "shal- 
low field"), while an additional 175 hours of observation 
were concentrated on the central portion of this field 
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Fig. 1. — The left and right panels show the sky coverage obtained during the January 2003 flight of Boomerang. The 
integration times per 3.4' pixel are shown for the shallow and deep data subsets described in the text. Galactic latitude contours 
of b = { — 30, —20} are shown for reference. Approximately 75 hours of observation are devoted to the shallow region, with an 
additional 120 hours concentrated in the central deep field. The cross-linking of the scans due to sky rotation is evident in the 
diagonal striping of the coverage, which has a « 20° opening angle. Known extragalactic point sources, shown masked in these 
plots, are excluded from the analysis. 



(the "deep region"). The deep field constitutes roughly 
~ 0.22% of the full sky. The exact distributions of in- 
tegration time for the data included in this analysis are 
shown in Figure Q] The B03 sky coverage is a subset of 
the region observed during Boomerang's 1998 flight. 

2.3. Inflight Calibration 

Both the calibration of the receivers and the charac- 
terization of the effective beam must be obtained from 
the inflight data. The former is due to the dependence of 
the bolometer responsivity on the radiative background, 
while the latter depends on the fidelity of the pointing 
reconstruction. 

The calibrations are obtained from the CMB itself; 
both the CMB dipole and degree-scale anisotropics pro- 
vide a well calibrated signal with precisely the desired 
spectrum. Scan synchronous noise limits the accu- 
ra cy of the d i pole c alibration to ~ 15%. As described 
in lMasi et alJ 12005 K cross- correlation with the WMAP 
data [Bennett et aL l|2003aD ] provides an absolute cali- 
bration uncertainty of 1.8%, including the uncertainty 
in the WMAP calibration. The B03 calibration is lim- 
ited by the relatively low signal-to-noise of the WMAP 
data at degree angular scales. The relative calibrations 
of the B03 detectors are obtained from cross-correlations 
of single channel temperature maps, and are accurate to 
0.8%. 

The error in the inflight pointing reconstruction is es- 
timated from the angular size of the extragalactic point 
sources in the B03 CMB fields. As indicated in Tabled 
the effective beams are consistent with the physical op- 
tics model of the optical system, convolved with a 2.4' 
r.m.s. Gaussian pointing jitter. Measurements of five 
bright sources in the shallow and deep fields provide an 
estimated 0.23' fwhm uncertainty on the width of the 
effective beams. 



3. ANALYSIS 

As a probe of the robustness of both the data and 
its analysis, the Boomerang team has implemented 
two fully independent analyses of the B03 time ordered 
data. The two pipelines, one centered in North Amer- 
ic a (NA) and one in Italy (IT), are descri b ed in detail 
inlMasi et alJ 120051) . IContaldi et all l|2005|U Jones et alJ 
l)2005lk lde"cTasDerisetHidll l)2005j) and references therein. 
In this section, we summarize the general approach to 
signal, noise, and angular power spectrum estimation. 

3.1. Signal/Noise Estimation 

Estimates of the signal on the sky and noise in 
the timestreams are obtained using an iterative pro- 
cedure si milar to that applied t o the B98 tempera - 
ture data iNetterfield et~atl l|2002fk iPrunet et all l|2001jk 
iRuhl et alJ <|2003ft ]. Least squares maps of the Stokes 
I, Q, and U parameters are generated from a combined 
analysis of the data from all eight 145 GHz PSBs, cor- 
responding to roughly 3 • 10 8 time domain samples. The 
data from each channel are weighted by their NET and 
combined during pixclization to decorrelate the linear 
Stokes parameters at each pixel. The Stokes parame- 
ter maps arc generated at ~ 3.4' resolution, resulting in 
approximately 377000 spatial pixels. 

The NA team divides the timestream into 200 noise- 
stationary subsets. The noise power spectra are calcu- 
lated from each of these chunks using the converged esti- 
mate of the signal derived from the full set of data. Each 
of these spectra are corrected for bias using an ensem- 
ble of signal and noise Monte Carlos. The IT analysis 
assumes stationarity of the noise over the course of the 
flight. While the NA approach is more general, the IT 
analysis benefits from a reduction in the sample variance 
of the noise estimate in the time domain. Each team em- 
ploys a self-consistent treatment of the correlated noise 
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Fig. 2. — In the top panel, we show a comparison of the 
temperature angular power spectrum derived from the B03 
data using the NA pipeline (xFaster, in red circles) and the 
IT pipeline, the latter using both uniform (green stars) and 
noise weighted (blue squares) masks. For the IT results, uni- 
form weighting of the combined shallow and deep data results 
in nearly optimal errorbars in the signal-dominated regime, 
while the l/\/N b 3 pixel weighting is required to achieve this 
sensitivity in the noise dominated regime. The NA spec- 
trum is derived using a hybrid auto- and cross- correlation 
technique which achieves nearly optimal errors across the full 
multipole range. In the lower panel, a compariso n of the an- 
gular power spectrum derived from the B98 data IRuhl et alJ 
(2003)], and the B03 results reported here. The 2003 cover- 
age is a sub-set of the region observed in 1998, so the signal is 
completely correlated between the two datasets. The binning 
of the B03 data in the lower panel is shifted by half of a bin 
relative to the top panel. 

between channels. Details regarding the tw o signal and 
noise estim ation pipelines are prov ided in Uones et alJ 
|2005|k and lde Gasperis et all l|2005j) . 

3.2. Power Spectrum Estimation 

Both the NA and IT analysis pipelines employ varia- 
tions of the pseudo-C^ Monte Carl o approach to powe r 
spectrum estimation introduced in iHivon et alJ l)2002|) . 
Each incorporates a polarized implementation of a gen- 
eral least squares signal estimator. The IT team uses 
a traditional pseudo-C^ approach using either \j \J Nobs 
or uniformly weighted masks to derive the tempera- 
ture power spectrum, while a cross-correlation tech nique 
similar to tha t desc ribed in iTristram et ah ( 2005|) and 
iPolenta et alJ (|2004|) is used for the polarization anal- 



ysis. Unlike these approaches, the NA estimator of the 
temperature and polarization power spectra is based on a 
hybrid auto- and cross-correlation technique adapted for 
multiple data sets characterized by partially overlapping 
sky coverage iContald i et alJ 1 20051)1. 

As described in IHivon et alJ l|2002f) . the Monte Carlo 
approach approximates the spherical harmonic transform 
of the heuristically weighted data, Cj T , as a convolution 
over the underlying power spectrum, Cj, T , with an ad- 
ditive noise term. Np. 



Cj T = " S ^2,Ku'Bf l FgiCj, T 



(3) 



The coupling kernel, Kw , is derived analytically from 
the transform of the weighted mask that has been ap- 
plied to the data. The window function, Bg, is the trans- 
form of the effective beam on the sky, and approximates 
the convolution over the Stokes parameters appearing in 
Equation^ The transfer function, Ft, is a symmetrized 
approximation of the combined effect of the scan strat- 
egy and time-domain processing on the sensitivity to a 
given angular scale. 

The transfer function is determined from an ensemble 
of signal-only Monte Carlo simulations, while the noise 
pseudo-power spectrum, Np, is obtained from an ensem- 
ble of noise-only simulations. Both the signal and noise 
simulations treat the flagged portions of the time stream 
in exactly the same manner as the actual data. 

3.3. The BOS Temperature Power Spectrum 

The Monte Carlo approach to power spectrum es- 
timation represents an approximate treatment of the 
noise covariance matrix of the map, and therefore gener- 
ally results in larg er uncertainties than do direct meth- 
ods |Borrilll l)1999|) ] . The sensitivity achieved by a Monte 
Carlo estimate of the power spectrum is determined by 
the properties of the noise and the weighting applied to 
the data. 

When analyzing a single map generated from the full 
set of B03 data, a complication arises in the selection of 
the weighting that is applied to the mask. In the limit 
of white (that is, pixel-uncorrelated) noise, the optimal 
power spectrum estimate of a noise-dominated sig nal is 
obtained by weighting the pixels according to 1/ y/Nobs- 
For the combined set of shallow and deep data, this 
scheme is therefore virtually equivalent to analyzing only 
the deep region. Conversely, in the signal dominated 
regime a uniform weighting will reduce (to the extent 
possible) the sample variance in the band power esti- 
mates. The choice of weighting reflects a trade-off be- 
tween the optimality of errorbars at low and high multi- 
poles. For the B03 data, the angular power spectrum is 
signal dominated at multipoles below £ < 900, as indi- 
cated in the fourth column of Table El 

The IT temperature power spectrum is derived from an 
analysis of the combined shallow and deep field data, us- 
ing an approach s i milar to th at applied to the B 98 data in 
INetterfield etaD J200l and IRuhl et all J200l . The IT 
analysis is carried out using either uniform or 1/ \J N b s 
pixel weighting. The agreement between the NA and 
IT analyses , as w ell as with the B98 results reported in 
IRuhl et all l|2003|) . is illustrated in Figure [3 These re- 
sults, representing two independent datasets and no less 
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Fig. 3. — The temperature angular power spectrum, CJ T , derived from the B03 data. The red points represent an alternate 
binning, and should not be interpreted as having statistical weight beyond that of the nominally binned data, which are shown 
with errorbars. The anti-correlations between neighboring bands range from the twelve to twenty percent level. The envelope of 
the beam uncertainty is indicated by ticks bracketing the one sigma errorbars for each of the bandpower estimates. As discussed 
in Section [3.51 these limits indicate the amount of "tilt" to the spectrum that is allowed by the beam uncertainty, and should 
not be interpreted as an additional uncorrelated error in each bin. The solid green line is the concordance ACDM model which 
best fits all published CMB data, including the B03 temperature and polarization results. The power spectra of the consistency 
checks described in the text are shown in the lower two panels. 



than three independent analysis pipelines, arc indicative 
of the robustness of the result. 

The NA pipeline addresses the problem of heuristic 
weighting of data with uneven coverage by treating the 
shallow and deep time ordered data separately. Each 
of these subsets, which have roughly uniform noise per 
pixel, exhibit partially correlated signal and statistically 
independent noise. A uniform weighting is applied to the 
shallow data, while the deep data are noise weighted. 

Estimates of the underlying power spectra are derived 
from the joint analysis of the two auto-spectra and the 
cross-power spectrum of the shallow and deep data, using 
a diagonal approximation to the quadratic Fisher ma- 
trix estimator. The inclusion of both the auto and cross 
power spectra allows nearly optimal errors to be obtained 
over the full range of angular scales, while computing the 
complete Fisher matrix in a self-consistent fashion for all 
the ban dpowers (temperature and polarization) simulta- 
neously |Contaldi et"aT1 lj2005j) ]. 

The full resolution B03 spectrum is displayed with the 
best fit conc ordance ACDM model i n the upper panel 
of Figure El MacTavish et~aTl ffiioftl }. The maximum 
likelihood bandpowers, the diagonal components of the 
correlation matrix, and the ratio of sample to noise vari- 
ance are presented in Table |5J As indicated in the fourth 
column of this table, the spectrum is sample variance 



limited at multipoles below £ < 900. Neighboring bins 
are anti-correlated at the ~ 12-20% level, the variation 
resulting from the exact form of the masks and the mul- 
tipole binning. 

The B03 data improve on the published measurements 
of the temperature angular power spectrum primarily 
over the third peak, at angular scales corresponding to 
the photon mean free path at the surface of last scatter- 
ing. All of the B03 temperature and polarization power 
spectra, inverse Fisher matrices, window functions and 
explanatory supplements are publicly available on the 
Boomerang web servers 1 . 

3.4. Internal Consistency Tests 

As a check on the internal consistency of the data, we 
perform two complete jackknife tests in which the TOD 
are divided in two halves, and /, Q, and U maps are gen- 
crated from each half independently. The power spectra 
are then computed from the difference of the resultant 
maps. 

The temporal jackknife divides the data into a first half 
and a second half, in which the shallow and deep scans 
are each divided in order to ensure roughly equal cover- 

1 http://cmb.phys . cwru. edu/boomerang/ index .html 
http : //oberon.romal . inf n. it/boomerang/b2k 
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Fig. 4. — The propagation of instrumental calibration errors 
to the bandpower estimates. The amplitude of the statis- 
tical uncertainties (sample variance plus instrumental noise) 
for each bin are shown with circles. The contributions from 
the relative calibration, s ± 0.8%, the the in-flight transfer 
function r ± 10%, the polarization angle, ip ± 2°, and polar- 
ization efficiency e ± 3% are shown for reference. The hor- 
izontal ticks indicate the effect of beam uncertainty. Note 
that each of these effects are highly correlated between bins, 
and therefore are not properly treated by quadrature addi- 
tion to the statistical uncertainty in a given bin. Since the 
correlation structure of these errors is known, they are more 
properly treated as nuisance parameters, and as such are 
mar ginalized over durin g the cosmological parameter estima- 
tion [Bri^^iOiLl J2002D 1. For the B03 parameter extraction, 
only the beam error and calibrat i on un certainty are treated 
in this manner |MacTavish et all (l200ll l. 



age of the s ky in e a ch sub set. As discussed in lMasi et alJ 
ll2005l and IJonesI J200oT) . B03 experienced a dramatic 
(~ 9km) loss of altitude over the course of the flight. 
The first-half /second- half test was chosen to provide a 
check which is maximally sensitive to systematic effects 
related to the altitude drop, such as responsivity drifts, 
a degradation in the accuracy of the pointing reconstruc- 
tion, or atmospheric contamination. 

The channel jackknife measures the difference between 
maps generated from the two halves of the focal plane. 
Each side of the focal plane accommodates two pairs of 
PSBs, allowing complete characterization of the three lin- 
ear Stokes parameters. The sky coverage, filtering, and 
data flagging of this jackknife test are nearly identical to 
that of the full data set. 

The noise properties, and therefore the inverse noise fil- 
tering of the signal, varies slightly from channel to chan- 
nel. When the data are processed in the time domain, 
this asymmetry, combined with the channel specific flag- 
ging of the time ordered data, results in a non-zero dif- 
ference spectrum even in a noiseless observation free of 
systematic effects. The temporal jackknife is similarly 
affected; the observations are not symmetric for the two 
subsets of data, resulting in subtle differences in the sig- 
nal processing. For both the channel-based and temporal 
divisions of the data, this offset is easily characterized in 
the Monte Carlo analysis by performing the same consis- 
tency test on an ensemble of simulations. The ensemble 
average of these simulated difference spectra are directly 
subtracted from those derived from the data. 

The high signal to noise ratio with which the temper- 
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The angular power spectrum of the B03 temperature 
data and the r esults of the consistency tests described 
in Section 13.41 The bandpowers, Cb £b(&b + 1)/27T and 
their errors are in units of [/iK 2 ]. The full set of temper- 
ature and polarization angular power spectra, and the 
associated covariance matrices and window functions, 
are publicly available on the Boomerang web servers, 
http : // cmb . phy s . cwru . edu/boomerang/ index . html 
and http : / / oberon . romal . inf n . it/boomerang/b2k. 
Neighboring bins are anticorrelated at the 12% to 20% 
level, depending on the multipole bin. The errors, ACb, 
are taken from the diagonal elements of the covariance 
matrix. In particular, the diagonal elements of the 
covariance matrix can be approximated by the formula 

Ac " - \j At ( c » + N ») > where Tb ^ 1 is 

a parameterization of the mode loss resulting from the 
filtering applied to the data. The Nb term represents the 
(binned and deconvolved) power spectrum of the noise, 
as determined from noise only simulations. Therefore, 
the ratio Cb/Nb indicates the relative contribution of 
sample and noise variance for a given bandpower. This 
quantity is shown for reference in the fourth column. 
The bandpowers obt ained from the two consistency 
tests described in 13.41 are given in the last two columns. 
We have subtracted the ensemble averaged power 
spectra derived from the jackknifes of the signal and 
noise Monte Carlos. The \ 2 to zero for the channel 
(temporal) jackknife is \ 2 = 20.8(55.5) for 24 degrees 
of freedom. 



ature anisotropics are detected is evident in the spectra 
of the differenced data, which are shown in the lower 
two panels of Figure The first of these is the spec- 
trum of the difference generated from the channel based 
division, the second is the spectrum derived from the 
temporal division. While the former (with a x 2 = 20.8 
for 24 degrees of freedom) passes the jackknife test, the 
latter clearly fails (% 2 = 55.5). The failure of the first- 
half/sccond-half consistency test, though significant, is 
at an amplitude which is small compared to the uncer- 
tainty from instrumental noise and sample variance in all 
but the last two bins, as shown in Table 
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Fig. 5.— A comparison of rece n tly pu blished CMB power spectra, Cj T , [the data shown are from iHinshaw et al.l (pf 
IKuo et all <I2004T ) : TReadhead et alj 1200411 . in addition to B03]. The data shown are derived from independent measurements 
that span two orders of magnitude in electromagnetic frequency (20-200 GHz), employ both single dish and interferometric 
instruments, and operate from terrestrial, balloon-borne, and orbital platforms. The remarkable degree of accord between such 
a diverse set of measurements is indicative of the maturity of the observational field. 



3.5. Propagation of Systematic Errors 

Uncertainties in the calibration and low level analysis 
of the data propagate in a complicated way from the 
TOD to the band power estimates. For the B03 analysis 
the principle uncertainties are the determination of the 
inflight transfer functions, the relative calibrations, the 
polarization efficiencies, and the effective beamsize. 

A very powerful feature of the Monte Carlo approach is 
the ease with which such effects can be modeled as they 
appear in the time domain, and properly accounted for in 
the analysis of the power spectra. In particular, we inves- 
tigate the impact of the uncertainty in the polarization 
leakage (e, ±3%), the PSB orientation in the focal plane 
(ip, ±2°), and the relative calibration (s, ±0.8%) that 
appear in Equations ^ and [21 Additionally, we probe the 
impact of the uncertainty on the inflight system transfer 
function (parameterized by a single-pole time-constant r, 
±10%) that is used to deconvolve each of the bolometer 
timestreams at the earliest stage of the analysis 2 . 

A simulated TOD is generated for each B03 channel 
using a realization of the best fit concordance ACDM cos- 
mology while randomly varying the instrumental param- 
eters about their nominal values. An ensemble of these 
TODs are processed through the NA analysis pipeline, 
and the distribution of the resultant bandpowcrs are used 

2 The bolometer transfer function is dependent on the inflight 
loading conditions and, for B03, is difficult to determine with the 



to propagate uncertainties on the instrumental parame- 
ters to systematic errorbars on the band power estimates. 
The results of this analysis, shown in Figure El indicate 
that B03 is not limited by systcmatics resulting from un- 
certainties in the characterization of the instrument. 

It is straightforward to analytically determine the im- 
pact of an uncertainty on the effective beamsize on the 
bandpower estimates. An estimate of the beam uncer- 
tainty is obtained from fits to the five brightest extra- 
galactic sources in the CMB field. In the limit that the 
pointing errors are Gaussian and isotropic, this misesti- 
mation will lead to a bias in spectrum equal to the ratio 
of the estimated and underlying window functions. This 
bias has a known spectrum 

-a\ (<5 2 +2«5) l(t+l) 



W'jW e 



(4) 



req uired accuracy fro m the inflight data itself. See I.Tonesj 
and lMasi et alJ I200fil) for further discussion of this topic. 



for an effective beam of width o~b = #fwhm/V8 In 2, and 
a fractional error, 5. The envelope of the two-sigma limits 
corresponding to a beam uncertainty of 8 ~ ±2% are 
shown as ticks bracketing the errorbars in Figure |3| 

The effects of both instrumental calibration and beam 
uncertainty are highly correlated from bin to bin, and 
have a known spectrum (Equation^, or are only weakly 
dependent on the underlying spectrum of the signal. As 
such, they are both properly treated as nuisance param- 
eters rather than independent contributions to the un- 
certainty in each bin. In the B03 parameter estimation, 
the calibration and beam uncertainty are tr eated in this 
manne r, following closely the prescription of lBridle et alJ 
(200J). 
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3.6. Foreground Contamination 

Galactic microwave emission can potentially contami- 
nate observations of the CMB. At 150 GHz, the thermal 
emission of interstellar dust grains is expe cted to be the 
domin a nt component of ga l actic e mission iBrandt et afl 
iITqqI : Ide Oliveira-Gostal lEffll : IMasi et alJ (120011) 1. 
Therefore any galactic contamination in the 143 GHz B03 
data is expected to exhibit some spatial correlation with 
existing infrared surveys. The spectrum of the galac- 
tic continuum emission is distinct from that of the CMB, 
allowing multifrcquency observations to discriminate be- 
tween the cosmological signal and local foregrounds. For 
this reason, B03 made simultaneous observations at 245 
and 345 GHz. 

Usin g the method described in lde Oliveira-Costa et al.1 
l)1999j) . we quantify the level of foreground contamina- 
tion in the B03 145 GHz data by cross-correlating the 
Stokes I parameter map with the Schlegel-Finkbeiner- 
Davis (SFD) dust map. The SFD map is a com- 
posite of_the_albsky DIR B E and IRAS lOOi i m sur - 
veys IHanser et alJ JTflflSh : INe 11g eha,uer et a,1J {Mi; 
iFinkbeiner et a,l.l ill 99QD : ISchlegel et a,1J l|1998D ]. 

We estimate the uncertainty on the correlation by 
applying the same analysis to fifty B03 noise realiza- 
tions. Before calculating the correlations, we resam- 
ple the DIRBE/IRAS data with the B03 scan pattern 
and filtering, and convolve both the DIRBE/IRAS data 
and the noise with a one-degree Gaussian. As a con- 
sequence, our results are insensitive to the zero point 
of the 100/xm maps. Using this procedure, we find no 
statistically significant correlation between the 145 GHz 
B03 and 100/xm data in either the B03 deep (b < —30) 
or shallow (b < —20) fields from which the CMB power 
spectrum is derived. 

The same analysis applied to the 245 GHz and 345 
GHz B03 data results in very clear correlations with the 
100 /im SFD map. The am plitude of these co rrelations, 
described in more detail in Ma si et aD ((2005) , provides 
an empirical scaling of the observed surface brightness of 
the 100 /xm dust to 145 GHz. When applied to the SFD 
map, this empirical scaling implies r.m.s fluctuations due 
to unpolarized dust of about 1 /xK in the B03 deep region 
at the angular scales of interest. This is roughly two 
orders of magnitude below the CMB contribution, and 
will result in a correspondingly negligible impact on the 
estimate of the CMB power spectrum. 

3.7. Features in the Temperature Power Spectrum 

A series of acoustic peaks is readily apparent in the 
B03 power spectr um. As has become common practice 
(see, f or example, Ide Bernardis et all l|2002|) ; IRuhl et aD 
OTOl : IPage et alJ II2003D 1. we calculate a model in- 
dependent characterization of the location, amplitude, 
and significance of the features in the power spectrum 
through a comparison of the goodness of fit of a parabola 
(or Gaussian) to each set of five contiguous bins 3 . In this 
work, we fit a Gaussian to the band powers in the vicinity 
of the first peak, and a parabola to the other features. 

3 While the result is relatively insensitive to the number of bins 
that are fit, as the subsets get larger than the characteristic size of 
the features in the spectrum, the reduced \ 2 clearly will degrade. 
Five bins is found to be the largest set that does not result in a 
poor goodness of fit statistic over the full range in I. 



As is well known, the lik elihoods of the band powe rs are 
not Gaussian distributed |Bond et all l|1998l I2000|) ]. We 
therefore transform to the "offset log-normal" variables 
whose likeliho od distributions are better approximate d 
by a Gaussian jBartlett et aD p000l) : lBond et all l|2000|) ] . 
The transformation is a simple one, 

Z h = \n{C b + Xb) 
where the noise offsets, Xb, are determined from the 
(binned and deconvolved) Ni as derived from the noise 
Monte Carlos. The inverse Fisher (covariance) matrix 
must similarly be transformed, 

F- 1 

(F z )- 1 = 

v >w {Cb + x b ){Cv + Xb>y 

In addition, the models to be fit arc binned into band- 
powers according to the same instrumental window func- 
tions, We, (shown in Figure EJ) which arc applied to the 
raw spectra 

where the logarithmic binning operator is defined as 

The parabolic model, 

Cf = C c (£-i ) 2 +C Q 

is similarly transformed using the offset log-normal ap- 
proximation, and the three dimensional likelihoods are 
calculated directly on a grid about each of the best 
fit locations. The A\ 2 contours for the curvature- 
marginalized likelihoods are shown in Figure 

In previously published B98 data, the significance 
of a detection has been determined by the curva- 
ture of the likelihood a t the peak of the distri- 
bution |de Bernardis et all {2002]k IRuhl et alJ (|200^ ]. 
However, as is evident in Figure the distributions are 
highly non-Gaussian. We therefore determine the sig- 
nificance of the detections from the amplitude of the 
marginalized likelihood at zero curvature. 

The first three peaks and three dips in the power spec- 
trum are detected with high confidence. Although the 
data favor a fourth peak in the vicinity of £ = 1055, with 
an amplitude of C1055 = 1020 [J.K 2 , the marginalized 
likelihood for the curvature parameter is not inconsistent 
with zero. For comparison, the same analysis was applied 
to the B98 data from the IRuhl et "aD (120031) release, as 
well as the binned first-year WMAP data [Bennett et alJ 
l|2003a|) ]. The results of all three analyses are compared 
in Table [21 and indicate a remarkable degree of consis- 
tency between the three independent experiments. 

The degree of concordance in temperature observations 
is further illustrated in Figure [5] which shows a com- 
pilation of the power spectrum estimates derived from 
four experiments, representing four very different exper- 
imental approaches, with observations probing nearly a 
decade in electromagnetic frequency. 

4. CONCLUSION 

We derive an estimate of the angular power spectrum 
of the CMB from the data obtained with the 145 GHz po- 
larization sensitive bolometers that flew on the January 
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2003 flight of Boomerang. The 245 and 345 GHz chan- 
nels place stringent limits on the level of foreground con- 
tamination. We characterize the systematic effects that 
result from various instrumental calibration uncertainties 
with Monte Carlo analyses, and verify the consistency of 
the data with jackknife tests. 

The B03 data are the first to be obtained using PSBs, 
which are identical in design to the polarized pixels in 
the Planck HFI focal plane. The sensitivity per resolu- 
tion element achieved in the B03 deep field is comparable 
to that anticipated by Planck at 143 GHz. The high sig- 
nal to noise of the temperature data results in a sample 
variance limited estimate of the power spectrum at mul- 
tipoles £ < 900. 

The B03 data presented in this work represent the most 
precise measurements to date of the angular power spec- 
trum between 600 < £ < 1200. In this regard, B03 plays 
a valuable role in bridging the gap between the all-sky 
WMAP survey at £ < 600, and the high angular resolu- 
tion data from terrestrial observations above £ > 1200. 
We characterize a series of features in the power spec- 
trum which extend to multipolcs £ > 1000, consistent 
with those expected from acoustic oscillations in the pri- 
mordial plasma in the context of standard cosmologies. 
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TABLE 3 

Features in the Temperature Power Spectrum 
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A comparison of the locations and amplitudes of the features in the temperature 
power spectrum derived from the B03, B98, and WMAP datasets. The values and (la) 
errors are obtained from the marginalized likelihood distributions directly. Values in 
parenthesis indicate a curvature parameter consistent with zero at 2a level, or greater 
(that is, a marginalized likelihood for the curvature parameter for which C(C C — 0) > 
2a). Note that this analysis has been performed on the binned WMAP data for 
better comparison with the Boomerang window functions. The full resolution WMAP 
spectrum provides constraints on the first peak and dip locations that are stronger than 
(but also consistent with) those reported her e. For the full resolution analysis of the 
first year WMAP data, see lPage etafl J2003h . 
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Fig. 6. — The left panel shows the CMB spectrum with the one and two-sigma Ax 2 contours for the fits shown in red (blue) 
for the features determined to have negative (positive) curvature. The likelihoods have been marginalized over the curvature 
parameter. The fourth "peak" only marginally favors negative curvature over a flat bandpower, and is not considered a detection. 
The right panel shows the window functions used in the generation of the band power estimates. 
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Fig. 7. — The marginalized likelihoods for the curvature parameter, C c , (at left) and the multipole, £o, (at right) of each feature 
in the temperature power spectrum. All three dips in the power spectrum are detected with high confidence, whereas only the 
first three peaks are detected with curvature significantly different than zero. The data around the first peak are fit with a 
Gaussian rather than a parabolic model, and the likelihood contours for this feature are not shown. 
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